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Abstract 

Using the Bethe ansatz, we obtain the exact solution of the master equa- 
tion for the totally asymmetric exclusion process on an infinite one-dimensional 
lattice. We derive explicit expressions for the conditional probabilities 
P(xi, . . . , xjv; t\yi, . . . , un] 0) of finding N particles on lattice sites x±, . . . , xn at 
time t with initial occupation y\, . . . , at time t = 0. 
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1 Introduction 



Driven lattice gases and particularly the one- dimensional asymmetric simple 
exclusion process (ASEP) have been intensively studied over the past decade 
for a variety of reasons []!]. The ASEP (for a definition see below) has been 
suggested already in 1968 as a model for the kinetics of biopolymerization 0. 
Two years later this process was introduced into the mathematical literature 
where it has received considerable attention in the context of interacting 
particle systems [|J. More recently the ASEP has been studied mainly by 
physicists as a model for polymers in random media and as a dynamical 
model for interface growth ||. It is also a discrete version of the noisy 
Burgers equation || and thus of interest for the study of shocks [0, || and for 
traffic models 11 . While there are exact solutions and a good understanding 
of the stationary behaviour of the system with (trivial) periodic and (non- 
trivial) open boundary conditions with injection and absorption of particles 
|T0| - []13| , exact results for the dynamics of the model are scarce (see e.g. 



Refs. @, [L4J] and, for more recent work, Refs. || y| |16| and references 
therein. It is the aim of this paper to present a new approach to this open 
problem by explicitly solving the master equation for the system defined on 
an infinite lattice. 

We study the ASEP with sequential updating. In the totally asymmetric 
version of this very simple model each lattice site can be occupied by at 
most one particle and particles hop with rate 1 to their right neighbouring 
site if this was empty. If it was occupied, the attempted move is rejected. 
This defines a Markov process with state space X = {0,1} 5 where a given 
configuration n G X is the set of occupation numbers Uk = 0, 1 with the site 
label k G S. Alternatively, if one restricts oneselves to studying the system 
with an arbitrary, but finite number of particles, one may define the process 
on Y = {{(/}}, {ki} , {ki, k 2 } } which is the collection of all finite subsets 
of S. In this case one identifies a configuration n which has particles on sites 
k\,...,kisr with the corresponding set in Y. Here we consider the system 
with finitely many particles N defined on an infinite one-dimensional lattice 
S = Z. 

A convenient presentation of the ASEP is in terms of a master equation for 
the probability P(Bn; t) of finding N particles on sites = {ki, . . . , /cat} G 
Y at time t. Defining the ASEP alternatively in terms of a master equation 
on the state space X rather than Y has turned out to be useful in previous 
work. In this case the stochastic time evolution is manifestly seen to be 
generated by the quantum Hamiltonian of a spin-1/2 Heisenberg chain. This 
suggests the use of the Bethe ansatz and the quantum group symmetry for 
the calculation of energy gaps (which give e.g. the dynamical exponent of 



the system) |Tj| [L8[] and certain time- dependent correlation functions |T7|, |19 



For the purposes of this paper using the coordinate representation Y of the 
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state space is more transparent, as in this case the master equation can 
be solved directly and explicitly using the coordinate Bethe ansatz |2"T| . 
While here we consider mainly the totally asymmetric case, the partially 
asymmetric exclusion process may be solved in the same way. We discuss 
this for the two-particle problem. 

In fact, even though the two-particle system may seem trivial, it exhibits 
already some of the characteristic behaviour of the driven system at finite 
particle density. In the undriven system (i.e. in the symmetric exclusion 
process) the diffusive spreading of a local perturbation does not depend on 
the overall density. In contrast, in the asymmetric case the center of mass of 
a local perturbation in an otherwise homogeneous finite background density 
is, from numerical work, known to spread super diffusively [[[]. This implies 
a divergent collective diffusion constant as the number of particles tends to 
infinity. There are, however, no quantitative results for a system with a small 
number of particles. As we will show here, in a system of two particles this 
behaviour appears as an increase of the (collective) diffusion constant which 
turns out to be the single-particle diffusion constant plus a term proportional 
to the square of the asymmetry in the hopping rates. 

The paper is organized as follows: Sec. 2 gives the main result of the 
paper. We formulate the master equation and present a (non-constructive) 
proof that the expression derived from the Bethe ansatz in Sec. 3 is indeed 
the solution. This is done because the proof is elementary and no reference 
to the Bethe ansatz is necessary. A reader not interested in the derivation 
of the solution may therefore skip Sec. 3 where the solution is constructed. 
As a simple application, Sections 4 and 5 focus on properties of the two- 
particle system. In Sec. 4 we briefly study the partially asymmetric process 
where particles are allowed to move both to the right and to the left, but 
with different rates. We restrict ourselves to the exact solution for the two- 
particle system, but also explain how to obtain the solution for the general 
iV-particle case. In Sec. 5 we obtain a very simple new result, which shows 
that already in the two-particle system the diffusive behaviour of the driven 
system is substantially different from the symmetric, undriven process. In 
Sec. 6 we present our conclusions. 

2 Solution of the master equation 

Let P(B N ; t) be the probability of finding N particles on the set of sites 
B N = {&!.,...,&#} at time t. When considering the probability P as a 
function of the coordinates ki we always assume this set to be ordered, ki < 
fcj+i V i. It is important to note that as a function of its arguments ki the 
function P is well-defined in 7, N , i.e. also for e.g. ki = ki + \ or ki > fcj+i. 
However in this domain P is not a probability. In other words, in the domain 



2 



&n = ki < k 2 < . . . < k^ C Z , the function P is the probability defined 
above, whereas in Z N \Qn it is defined by the master equation below, but 
is not a probability 

For the totally asymmetric exclusion process as described in the intro- 
duction, P(E>n; t) defined on Z N x [0, oo) satisfies the master equation 

j t P(ki, ...,k N ;t) = P(h-1, ...,k N - t)+. . .+P(k x , . . . , k N -l; t)-NP{k u ...,k N ;t). 

(2.1) 

This has to be supplemented by boundary conditions in Z . If any two 
neighbouring arguments ki, k i+1 are equal, P satisfies 

P(h, ...,h, k i+ i = h,..., k N ; t) = P(h, ...,k h k i+1 = ki+1, . . . , k N ; t) V t > 0. 

(2.2) 

This boundary condition expresses the exclusion interaction. This is easy to 
see in the simplest case of two particles. Then ( |2~T| ) and (|2.2|) defined on Z 2 
read 

j t P{K k 2 ] t) = P(h - 1, k 2] t) + P(h, k 2 - 1; t) - 2P(h, k 2 ; t) (2.3) 

and 

P(k,k;t) = P(k,k + l;t) V k and t > (2.4) 
which is equivalent to the following equations with P restricted to fl 2 

^-P{kx, k 2 - 1) = P(ki - 1, k 2 ; t) + P{k t , k 2 -l;t)- 2P{k x , k 2 ; t) ]lk 2 -h>l 
at 

(2.5) 

and 

j t P(h, k 2 ; t) = P(h - 1, k 2 ; t) - P(k u k 2 ; t) if k 2 -k 1 = 1. (2.6) 

The second equation (|2.6j ) expresses that due to exclusion the configuration 
(k, k+1) can be reached in a single step only from the configuration k— 1, k+1 
and be left only in a single way (which is by moving to (k, k + 2)). Extending 
the range of validity of (|2.5|) to all Z 2 requires adding P(k\, k 2 — 1; t) — 
P(ki,k 2 ',t) to (|2.6|). However, by demanding that ( |2.4|) holds for all times, 
this is equivalent to adding 0. Thus (|2.5|) and ( p.6|) remain unchanged, i.e., 
the two sets of equations have the same solutions in the "physical" domain 
Q 2 . It may seem more natural to use the second formulation of the master 
equation which is in the iV-particle case the restriction of the validity of 
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the master equation ( |2.1| ) to the domain Qn- One could then obtain a set 
of equations equivalent to (|2.1|) , (|2.2|) by replacing the boundary condition 
( j2.2| ) by including appropriately chosen Kronecker delta-functions in ( |2.l|) . 
However, it turns out that solving the equation is more straightforward and 
transparent in the formulation ( |2.1| ), (|2.2j). 

We finally note that with specified initial condition An = {h, . . . ,In} £ 
Y, i.e., 

P(B n ;0)=6 An , Bn (2.7) 

the probability P(B N ; t) becomes the conditional probability P(B N ; t\A N ; 0) 
and thus a complete solution of the problem. Also as a function of the 
arguments l{ the function P is a probability only in the domain i x < . . . < 
Let us now introduce the function 



k=0 



where the binomial coefficient and the factorial are defined by the T-function, 
i.e. a! = T(a + 1) and 

T(a+1) . . 

(2.9) 



V r(6 + i)r(a-& + i)' 



In what follows we shall need only p,n G Z and t G [0, oo). We list some of 
the properties of F p (n; t): 



(1) For integer p < 0, F p reduces to a finite sum 

bl 



In particular, 

F (n-t) = -e- t (2.11) 
n\ 



(2) For the time derivative one finds 



jFM t) = F P-i( n - 1; *) = - 1; <) - ^ P (^; 0, (2.12) 
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(3) and for the time integral one gets 

dtF p (n-t) = F p+1 (n+l;t)-(- n - p 1+p ) = £ F p (k; t) - (2.13) 

k=n+l 



(4) At time t = one has 

limF p (n;t)=(-; + _7 1 ) (2.14) 

which vanishes for n > 0. Now we state the main result of this paper: 

Theorem: Let F(Bn, An] t) be the N x N matrix with matrix elements 
Fij = Fi_j(ki - Ij; t). Then 

P(B N ; t\A N ; 0) = det F(B N , A N ; t) (2.15) 

is the solution of the master equation ( |2.1|) with boundary condition (|2.2|) 
and with initial condition (127 



The theorem states that the conditional probability of finding iV particles 
at time t on B^ C fijy if initially (at time t = 0) they had been on sites 
An C Qn is given by the determinant Q2.15| ). How this result was derived 
is explained in the next section. Here we give a proof of the theorem which 
is independent of this construction. First we show that the determinant 
is a solution to the master equation ( ^.1[ ). Then we show that it satisfies 
the boundary condition ( j2.2|) and finally we prove that it satisfies the initial 
condition (|2.7|). 

Proof: (i) Because of the factor e~ l in the functions F p (n;t) the matrix F 
may be written e~ l F and one gets 

4-detF = e- m ^-detF - N det F. (2.16) 
alt dt K ' 

This accounts for the term —NP on the r.h.s. of (]2.1|). The time derivative 
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of the determinant of F may be written 



d 
dt 
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NN 



The matrix elements in row i are -^-2(^—^2), ■ ■ ■ , Fi^ N {ki—l N )) 

and their time derivatives are — 1 — 1%), Fj_ 2 (fcj — 1 — 1 2 ), ■ ■ ■ , Fi_ N {ki — 

1 — In)~)i see ( P-12|) . Thus each determinant in the time derivative of F 
contributes exactly one of the terms P(k\, . . . , ki — 1, . . . , fcjv; t) on the r.h.s. 
of ( |2.1| ), i.e. ( |2.15| ) satisfies the master equation ([2.11). 

(ii) In order to show that ( |2.15|) satisfies the boundary condition (|2.2| ) we 
note that according to fl2.12| ) each column in F may be written 



[Fi-i(ki - h), Fi_2(ki -l 2 ),..., Fi_ N (ki - l N )) 
= (F^ih + 1-/1), F^ 2 {ki + 1 - f 2 ), ... , F^ N (h + l-l N )) 
+ (Fi(ki + F^ki + l-k),..., F i+1 _ N {ki + l-l N )) (2.17) 
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Assume now that ki+i — ki + 1. This gives for the r.h.s. of ( |2.2| ) 





F {kx-h;t) 


F-iih - 


kit) 


. F X -. N {ki - l N ;t) 






Fi-^hi - fat) 


F o(k- - 




P. M (U. _ /. r - f) 






Fi(k i+1 -h;t) 


Fi_i(A; i+ i 


-hit) . 


■ F i+ i_N(k i+ i — 1^ 


t) 




F N „i(k N - hit) 


F N _ 2 (k N 


-hit) . 


■ F (k N -l N ;t) 






F-iih-h 


It) 


■ F 1 _ N (ki - l N ; t) 






^(ki + l-hit) 


F^ 2 {h + 1 


-hit) . 


■ Fi_N{ki + 1 — In] 


t) 


Fi(k i+1 -h;t) 


Fi_i(k i+1 - 


hit) . 


• ^i+l-Ar(&i+l — In 


t) 


F N ~i{k N -h;t) 


F N _ 2 (k N - 


hit) . 


■ F (k N -l N ;t) 





+ 



Foih-hlt) F-xih-hlt) 

Fiih + l-hit) F^tih + l-hlt) 

Fi(k i+1 - hit) Fj_i(A; i+ i - h] t) 

F N -i{k N -h;t) F N _ 2 (k N -h]t) 
Foih-hlt) F^ih-hlt) 



F^ih + l-hit) F^ih + l-hit) . 
Fi(k i+1 - hi t) Fi_i(A; i+1 - h] t) 

F N _ x {k N - hit) F N _ 2 (k N ~ hit) 



.. Fi_jv(A;i - l N ;t) 

.. F i+1 _ N (ki + 1 - l N ;t) 
■ ■ F i+ i_ N (k i+1 — l N ; t) 



F (k N -l N ]t) 

■ F^ N (ki - l N ] t) 

. Fi_ N (ki + l-l N ;t)^ 

■ F i+ i_ N (k i+ i — l N ; t) 

■ F (k N -l N ;t) 



.18) 



The second determinant in the sum on the r.h.s. of this equation vanishes 
since two rows are identical. It remains the first determinant in the sum 
which is equal to the l.h.s. of (|2.2|). 

(iii) It remains to show that (|2.15|) satisfies the correct initial condition 
Q2.7|) . We first assume that k\ > h. This implies ki > h since k± < . . . < k^. 
At t = 0, all the matrix elements Fn = Fi_\(ki — h;0) vanish, see ( |2.14| ), 
and hence the determinant vanishes. Now we assume k\ = h- Then in the 
first column F n = 1, but the other matrix elements F a in the first column 
are still zero, since all ki are larger than h except k\. The determinant of F 
is therefore equal to the determinant of the matrix F^ obtained from F by 
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omitting the first row and the first column. 

In the next step assume now first that k 2 > l 2 . Repeating the previous 
step gives det F(Bn, An] 0) = S^^S^,^ det F^ 2 \ Iterating this procedure N 
times finally gives det F(Bn, An] 0) = 5a n ,b n - D 

An integral representation of the terms in the determinant is derived in 
the next section (|3.14j) . 



3 Bethe ansatz solution 

For a derivation of the solution the master equation (|2.l|), (|2.2|), (|2.7|), we 
first turn the differential equations (|2.1| ) into an eigenvalue problem by the 
ansatz P(Bjy;t) = e~ et P(BN)- In order to solve for the resulting difference 
equation we follow the strategy employed by Bethe for the solution of the 
isotropic Heisenberg spin chain p(| and extended by Yang and Yang |21] to 



the anisotropic spin chain. Rather than using hi, U for the integer coordinates 
of the particles we shall use in this section the notation Xi, jji G Z. For 
momentum labels we shall use pj. 

First we consider N — 1. The resulting equation is 

eP(x) = -P(x - 1) + P(x) (3.1) 

which is readily solved by P(x) = e ipx with p e [0, 2tt). This gives for the 
"energy" 

e p = 1 - e~ ip (3.2) 

and P(x;t) = J^ n dpe~ 6pt f(p)e ipx . The initial condition P(x; 0) = 5 XtV deter- 
mines f(p) = e~ ipy /(27r) and finally yields 

P(x;t\y;0) = — [ dpe-^- ipv e ipx (3.3) 
27r Jo 

fx-y 

e- f (3.4) 



(x-y)\ 

F (x-y;t) (3.5) 



For N = 2 one has to solve 



eP(x 1 ,x 2 ) = -P( Xl - l,x 2 ) -P(x 1 ,x 2 - 1) + 2P{ Xl ,x 2 ) (3.6) 
P(x,x) = P{x,x + 1) Vx (3.7) 



S 



The first equation is solved by Bethe's ansatz 



P(x u x 2 ) = A 12 e ipiXl+ip2X2 + A 21 e ip2Xl+ipiX2 (3.8) 
with arbitrary constants Aij(pi,p 2 ) and gives 

€ Pl,P2 = e pi + € P2- (3-9) 



The second equation (|3.7|) fixes the ratio Si 2 = Ayij A 2X . Inserting flcTq) gives 



s * = < 3 - 10) 

The range of values pi and p 2 may take needs some discussion. In 
the usual Heisenberg quantum chain one finds a bound state in the two- 
particle sector, i.e. a state with complex momenta p\ )2 = u ± iv. This 
is a solution for vanishing v4 12 or vanishing A 2 i in which case the wave 
function decays exponentially in the distance x 2 — x\ (see next Section). 
Here there is no non-zero p for which either A\ 2 or A 2 \ vanish and hence 
no bound state. We conclude that p\,p 2 G [0, 2tt) and P(xx,x 2 ;t) = 
J dpxj dp 2 e- {t ^ + ^) t f( y p u p 2 )( e ipiXl+ip2X2 + S 21 e ip2Xl+ipiX2 ) is the general so- 
lution of ( p.l| ) with boundary condition (|2.2|). (For obvious reasons we define 
S 2 i = S 12 .) 

In order to satisfy the initial condition Q2.7] ) one has to determine f(pi,p 2 ) 
and discuss the pole resulting from the integration over S 2 i- Assuming that 
the particles were initially at sites y±, y 2 it turns out that choosing f(pi,p 2 ) = 
e -ipiyi-iP2V2 anc [ defining the position of the pole in S 2 \ by p\ — > p\ + iO gives 
the correct initial condition P(xi,x 2 ;0) = S xl}yi S X2ty2 . This gives 

P{x u x 2 \t\yi,V2\G) = / dpi / dpae-^+^^Hww-iww 

(2vr) 2 J J 

/ 1 _ J-P2 . \ 

\ l^vl e ) [6AL) 
= F (xi - yi;t)F {x 2 - y 2 ;t) 

-F-i(x 1 -y 2 ;t)F l (x 2 -y 1 ;t) (3.12) 

= det F(B 2 ,A 2 ;t) (3.13) 
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with £> 2 = {xi, X2} C Q2, A 2 = {yi, JJ2} C Q2 and the position of the pole in 
( [3. 11| ) defined as discussed above. 



In the same way one continues for N > 3. One constructs a superpo- 
sition = ^2 A il; _ tiN exp (ip^Xi + . . . + ipi N xw) of iV-particle plane waves 
with all iV! possible permutations of the wave numbers pi and with ampli- 
tudes Ai ly __ } i N . The ratio of any two amplitudes for plane waves where two 
momenta Pi,Pj are interchanged is Sij, in the same way as in the two-particle 
case (|3.10|) . This takes care of the boundary condition fl2.2|) when (any) two 
particles are on nearest neighbour sites. The crucial point is that for higher 
number of particles there are no new constraints from the boundary condition 
when more than two particles are on adjacent sites. This can be seen by not- 
ing that satisfying the boundary condition for any given pair is independent 
of the coordinates of the remaining particles. So one constructs the Bethe 
wave function by starting from e i Pi x i+-+ i PN^N w ith amplitude A12...N = 1 
and then performing all possible permutations of the momenta. For each 
permutation — > one multiplies with a factor Sji as in (|3.11[ ). The 
total "energy" e corresponding to such a wave function is the sum of the 
single particle energies e pli ..., Pjv = J2iLi e Pi- The initial condition (|2.7f ) deter- 
mines the overall normalization of the wave function and the position of the 
poles arising from the integration over the various S^ appearing in the wave 
function. It is satisfied by the choice f(pi, . . . ,Pn) = e~^ piyi+ - ipNVN ^ and by 
placing the poles in the same way as in the two-particle case, i.e., by setting 
Pi — > Pi + iO in the denominators. 

This construction provides an integral representation of the N\ terms 
appearing in the determinant ( 2.15|) . Therefore the solution of the master 
equation may be written 



N ^ ^2tt 

P{x x , ...,x N ] t\x u ...,x N ;0) = Y[— dp,' x 



= n s y o "n 

ty pi ,..., PN (x 1} ... } x N ) (3.14) 
with the Bethe wave function \1/ as defined above. For three particles it reads 

\T/ ( rr <r rr \ — „ipiXl+ip2X2+ip3X3 i Q J,p2X\+ip\X2+ip3X3 

V Pl,P2,P3V X l> X 2; X 3j — C T ^21 e 

I C C e iP2Xl+ip3X2+ipiX3 _|_ C C „ipzxi+ip2X2+ipiX3 
' 21 31 ' 21 31 32 

_|_ Q^g^ e ip3Xi+ipiX2+ip2X3 _|_ g^ e ipixi+ip3X2+ip2X3 
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4 The partially asymmetric process 

Using the Bethe ansatz one may also solve for the partially asymmetric pro- 
cess where particles are allowed to move with rate Dl to left and with rate 
D R to the right. The strategy is the same as for the totally asymmetric case 
discussed in the preceding section. The main difference in the analysis is the 
occurrence of bound states in addition to the continuum. We shall discuss 
in some detail only the one- and two-particle systems. This is, in principle, 
sufficient to construct the general ./V-particle solution. 

The case of a single particle can be copied with little modification from 
the previous section. It is convenient to introduce the asymmetry q = e a and 
the time scale D by 

q = yir L 

a = In q 
D = ^D R D L 

After separating the time dependence the master equation reads 

eP(x) = -D R P(x - 1) - D L P(x + 1) + (D L + D R )P(x) (4.4) 

which is readily solved by P{x) = e ipx with p e [0, 2ir). This gives for the 
"energy" 

e p = D R (l-e- ip ) + D L (l-e ip ) (4.5) 

and P(x;t) = dpe~ ep f(p)e tpx . The initial condition P(x; 0) = 8 x>y deter- 
mines f(p) = e~ tpy /(27r) and finally yields 

P(x;t\y;0) = — f" dpe'^-™ e ipx (4.6) 
2tt J 

= e- {q+g ~ 1)Dt q x - y I x _ y {2Dt) (4.7) 

where I n (r) is the modified Bessel function. The representation of (|4.6|) 
in terms of the Bessel function ( |4.7|) is obtained by an elementary contour 
integration. It is easy to verify that both expressions (|4.6|) and ( |4.7| ) sat- 
isfy the same differential-difference equation with the same initial condition 
P(x;0\y;0) = 5 x , y . 
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(4.1) 

(4.2) 
(4.3) 



Note that the Bessel function diverges asymptotically ~ e 2Dt which is not 
sufficiently fast to cancel the prefactor e~^ +9 ^ Dt in (|4.7| ). The interpretation 
of this observation is that the probability of finding the particle at site x with 
x kept fixed decays exponentially with an inverse correlation time or "energy 
gap" Ae = D(q + g" 1 — 2). However, by going into a comoving frame with 
velocity v = D^q—q^ 1 ) = Dr—Dl (see Sec. 5), i.e. by studying the behaviour 
of the distribution around x' = x + vt, one finds the usual algebraic, diffusive 
behaviour. 

For N = 2 one has to solve 

eP(x 1 ,x 2 ) = -D R (P(x 1 -l,x 2 )+P{x 1 ,x 2 -l)-2P(x 1 ,x 2 )) 

-D L (P( Xl + 1, x 2 ) + P( Xl , x 2 + 1) - 2P( Xl , x 2 )) (4.8) 
p , , ^ DrP(x,x)+D l P(x + 1,x + 1) 

P(x,x + 1) = — — Vi (4.9) 

The first equation is solved by Bethe's ansatz 

P(x u x 2 ) = A 12 e ipiXl+ip2X2 + A 21 e ip2Xl+ipiX2 (4.10) 
with arbitrary constants Aij(px,p 2 ) and gives for the two-particle energy 

e Pl,P2 = € Pl + € P2- (4-11) 



The second equation ( |4.9|) fixes the ratio 5i2 = A± 2 /A 2 i. Inserting (|4 . 1 0|) 
gives 

s Dr + D L e ipi+ip2 - (Dr + D L Y ipi f v 

12 ~ D R + Die ipi+ip ' 2 - (D R + D L )eF» 1 ' ' 

which depends only on the momenta pi j2 and the asymmetry q. 

As discussed in the preceding section here one finds besides the continuum 
Pi)f>2 £ [0, 27r) a solution corresponding to a bound state. To see this, we set 
Pi = u — i(a — v) , p 2 = u — i(a + v) with u, v real and a = \nq. Clearly, in 
order to obtain a wave function which decays exponentially in x 2 — X\ > 0, 
either A% 2 or A 2 i must vanish. Choosing 

e- = < 1 (4.13) 

cosh a 

gives A l2 = D R + D L e ipi+ip2 - (D R + D L )e ipi = 0. The wave function QOTD 
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reduces then to a single expression 



P(x u X 2 ) OC q ^+^ e iu(x 1 +x 2 ) e -v(x 2 -x 1 )_ ( 414 ) 

Using ( 4.11| ) one obtains for the "energy" of this state 



cos 2 u 

e u = 2D cosh all =- ) . (4.15) 

cosh a 



Thus the bound state has a non- vanishing energy gap Ae = 2D sinh 2 a/ cosh a 
for any asymmetry a ^ 0. 

For the determination of the solution of the master equation with initial 
condition S xim 8 X2iV2 one proceeds in a way analogous to the totally asym- 
metric case. First we note that choosing A 12 = 1, multiplying the time- 
dependent wave function by exp (—ipiyi — ip 2 y 2 ) an d integrating p\ and p 2 
from to 2tt gives, at time t = the correct initial value S xltyi S X2t y 2 plus a 
non- vanishing term from the reflected wave proportional to S 2 \. This term 
needs to be cancelled by an appropriate choice of the amplitude of the bound 
state contribution. This may be determined by using center of mass coordi- 
nates R = x\ + X2, Rq = yi + yi and relative coordinates r = x 2 — x\ > 0, 
r o = V2 — Vi > 0. Solving first for the .R-dependence of the master equation 
(which is trivial) one obtains then a lattice diffusion equation in a single co- 
ordinate r with partially absorbing boundary condition. This equation was 



solved in [22| and one finally finds that 

P(x 1 ,x 2 ;t\y 1 ,y 2 ;0) = d Pl £ tfoe-^i^Mwi/iHwia x 

^ e ipiXi+ip2X2 _|_ ^^ e ip2Xi+ipix 2 j 

+- [* du (e 2v - 1) q R - R ° e iu{R - Ro) e- v{r+ro) et446) 
n Jo 

— pcont _|_ pbound ^ 17) 



solves the master equation of the partially asymmetric process with two par- 
ticles and initial condition P(xi, x 2 ; t\yi, y 2 ; 0) = S Xl>yi S X2>y2 . The first piece 
pcont - n ^ e sum is the contribution from the continuum of states (|4.10 ), 
whereas in the second piece p bound one recognizes the contribution from the 
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bound state. Using the identity 



1 f 2n 

fi- n I n {2r) = -f- / dpe ipn+ ^ eW+ ^ le ~' P)T (4.18) 

27T ./n 



and expanding the denominator of S21 in a geometric series in z(pi,p2) = 
(Df i e~ ip2 + Di J e' ipi )/(Dji + D^) one may rewrite ([OB) in terms of a sum 



of products of two modified Bessel functions. With (|4.18|) the bound state 
contribution takes the form 

pbound = 1 ^ due iu(R-R ) A _ £2\ £.r+r ~2 e ~(l-e)(D R +D L )t (4.19) 
TT Jo 

where £(w) = z(u, —u). 

The general iV-particle problem is solved by the Bethe ansatz for N par- 
ticles and by determining the various bound state contributions. A way of 
determining the contributions from the bound states in the general case is by 
an appropriate contour integration in the complex ki planes which includes 
the poles of the reflection coefficients Sij. These poles give rise to the bound 
state contributions in a general iV-particle problem. 



5 Diffusion of two particles 

Here we want to study how the exclusion interaction affects the diffusion of 
two particles. In order to get an understanding of what is happening we 
put the particles at time t = on lattice sites y\ — — 1 and %)2 = 1 and we 
study the moments of the density distribution at time t. This describes the 
diffusive broadening of an initially spatially concentrated density. 

We introduce expectation value (n x ) which is the probability of finding 
a particle on site x at time t. The moments of this density distribution 
may be obtained from the Fourier transform p(q) = ^2 x e~ iqx (n x ) by taking 
derivatives w.r.t q. Here we are interested in 

N = p(0) = 2 (5.1) 

w = ^p'(o) = \J2 X ^ ( 5 - 2 ) 

X 

(X 2 ) = -^p"(0) = ^J> 2 K> (5.3) 
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from which we shall calculate the asymptotic drift velocity v and the asymp- 
totic collective diffusion constant A defined by 



v = hm^(X) (5.4) 

t^oo at 

A = lim ± ((X 2 ) - (X) 2 ) (5.5) 
t^oo at 

It may be useful to remind the reader of what these quantities are in 
case of non-interacting particles. This allows for a comparison of the in- 
teracting and the non-interacting system. For a single particle or for two 
non-interacting particles the density satisfies the diffusion equation 

^(n x ) = D R {n x ^) + D L (n x+1 ) - (D R + D L )(n x ). (5.6) 
at 

Integrating this equation gives after a short calculation 

v = D R -D L (5.7) 
A = D R + D L (5.8) 

which is trivial in the sense that this is just a way of defining the driven 
non-interacting process. The non-trivial point is the determination of these 
quantities for the system with exclusion interaction. 

First we note that for the ASEP the density satisfies the continuity equa- 
tion 

j t (n x ) = Ux-i) - (j x ) (5.9) 

with the current 



(jx) = D R {n x (l - n x+1 )) - D L ((l - n x )n x+1 ) (5.10) 

Thus one gets v = J2 x (jx)/2 = ( D R ~ D L)J2 x {( n ^) ~ (n x n x+1 ))/2 = 
(D R - D L )(1 - ExK^+i)/ 2 )- Using (fO|) one may write J2 x ( n xn>x+i) = 
J2 x P(x, x + l;t) = E x (D R P(x, x; t) + D L P(x + 1, x + 1; t))/(D R + D L ) = 
^2 x P(x,x;t). Thus with the initial condition considered above one gets 

(X) = (D R -D L )(^t-~^dTj2 Pfa *) r\ - 1, 1; 0) J . (5.11) 
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Using the continuity equation the diffusion constant may be written 
A = E*(z + 1/2) (jx) - 2v(X) = v + 2D L + 2(D R - D L - v)(X) - (D R - 
D L ) J2 x x(n x n x+1 ). Now (|5~TT|) and Q lead to 



A = lim {D R + D L + 



t^oo 



(D R -D L )*Y,P(x,x',t\-l,l;0) 



1 f* 



~(D R - D L ) 



Et^-^^y )P(x,x;t|-l,l;0) 



(5.12) 



Therefore in order to determine t> and P and has to calculate Pq = 
P(x, x; t\ — 1, 1; 0) and Pi = xP(x, x; t\ — 1, 1; 0). Each of these quan- 
tities can be split into three contributions arising from the two contributions 
from the continuum and the bound state contribution in (|4.16 ). 

There is no contribution from the bound state to 52 P(x, x; t\ — 1, 1; 0) 
and defining r = 2(D R + D^)t one finds 



P = e- T (/ 2 (r) + / 1 (r)) 



(5.13) 



where the term proportional to I\ is the term arising from the reflected wave 
proportional to S 2 i- 

The calculation of Pi is slightly more involved, but still straightforward, 
and gives 



Pi 



2(D R + D L 



-e- T {D R h-D L h)+e- T I 2 



T 



2(D R + D L 
Dr-Dl 

D R + Dr 



-e- T (D R I 2 -D L I ) 



Dr-Dl 
2(D R + D L ) 



D 



R 



D R + D L 



e- T h 



(5.14) 



where the last piece in the sum comes from the bound state and the argu- 
ments of the Bessel functions are all r. 

Putting everything together and taking the limit t — > oo in ( |5.12|) finally 
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yields 



D R - D L (5.15) 

Dr + Dl+ (J^zJ^(1-1\. (5 , 6) 



D R + D L \2 7T 

On comparison with the results ( |5.7| ), ( |5.8| ) for non-interacting particles one 
notices that the exclusion interaction alone does not change the collective 
two-particle diffusion constant. In the undriven system one has A = Dr+Dl 
as in the non- interacting system. In the presence of the drift, however, A 



increases to the value (|5.16|) . This result was obtained for particles placed 
initially at sites = —1, 1. It is however valid for any (finite) initial 
separation tq. After a time to 3> r o the details of the initial condition are 
washed out. 



6 Conclusions 

The main result of this paper is the solution ( |2.15| ) of the master equation for 
the asymmetric simple exclusion process. This solution allows for a complete 
description of a system of finitely many particles. As a simple example we 
have investigated the collective diffusion of two single particles. We found 
that the diffusive broadening of the density profile in the driven system ( |5.16| ) 
is faster than both in the undriven and in the non-interacting case. 

The solution of the master equation may also be used for the analysis 
of quantities in systems with finite density. Consider e.g. a system with 
constant non-zero density p with a local inhomogeneity such as a lattice site 
y where initially the density is (n y ) = 1. In such a situation it would be 
interesting to study the time evolution of the density profile (which gives the 
dynamical structure function) or the temporal behaviour of density correla- 
tions. Using the exact solution one can obtain an exact expansion of these 
quantities in powers of p where the n th power is obtained by solving the n- 
particle problem. This can be seen as follows: Suppose one wants to calculate 
the time-dependent density profile p x (t) = (n x (t)) up to second order in the 
background density p. The time derivative of the two-point correlation func- 
tion involves a three-point correlator which is of order p 3 and which therefore 
may be neglected in the desired second order approximation. Leaving the 
three-point correlator out results in a differential-difference equation for the 
two-point correlator which is identical to the two-particle master equation 
( |2.1|) with boundary condition ( |2.2|) . Thus one can calculate (n x n y ) and 



then by summing up two-point correlators one gets (n x ) up to order p 2 . For 
a third order approximation one considers the three-point correlation func- 
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tion. If one neglects fourth order correlators, it satisfies the three-particle 
master equation. Summing up three-point correlators yields (n x ) up to order 

A by-product of Sec. 5 is the explicit solution of the two-particle problem 
for the symmetric exclusion process. It is interesting to recall that in the 
symmetric case D R = D L the conditional probability P(xi,x 2 ;t\yi, 3/2! 0) 
determines not only the behaviour of the two-particle system, but also the 
behaviour of various time-dependent density-density correlation functions in 
iV-particle systems, viz. the equal-time two-point correlator (n Xl (t)n X2 (t)) 
for an arbitrary initial state ||, the two-time correlator (n Xl (ti)n X2 (t2)) for 
an arbitrary initial state |23[] and the (time-translationally invariant) four- 
point correlator {n Xl {ti)n X2 {t2)n X3 {tz)n Xi {ti)) averaged over the stationary 
distribution [23]. The contribution of the bound state to pair diffusion was 



discussed in [24]. It would be interesting to investigate in detail how the 



bound state effects the behaviour of these correlators in the various space- 
time regimes of the symmetric diffusion process. 

One puzzling problem is the behaviour of the ASEP in the presence of 
a blockage [^, E71, i.e., a bond in the lattice where particles hop with 
rate r 7^ 1. Numerical and analytical studies seem to indicate that the 
steady state current in a finite, half-filled system approaches its maximal 
value j max = 1/4 already at a surprisingly small defect hopping rate r ^ 0.8 



25 , p9| . This raises the question of a non-analyticity in the current j(r) 
for r < 1 which only an exact calculation of the steady state current can 
resolve. A perturbative expansion in r has been performed up to sixth order 



using computer algebra |28j. Each coefficient in the expansion is a rational 
number with numerators and denominators rapidly increasing with the or- 
der. Many non-trivial exact steady state properties of exclusion processes 
have been obtained by exact calculation for e.g. small system sizes, then 
guessing the general structure, and finally proving that the exact expressions 



obtained in this way are correct |1G, 11, 30, p]], B3, 33]. Unfortunately, ap 



plying this strategy to the perturbative coefficients obtained in seems 
hopeless. However, a systematic perturbative expansion of the current may 
be performed using the solution (|2.15|) of the master equation where the n th 
order in r is obtained by the solution of the n-particle problem. This gives 
the coefficients of r n as a sum of n! fractions with slowly increasing denom- 



inators and numerators rather than the single fractions given in Ref. |]28 
Therefore there is some hope that one might be able to guess a pattern in 
this sequence of fractions. Clearly this is a somewhat speculative suggestion 
which very well may turn out to be useless. But it seems worth trying to 
obtain exact results in this way. 
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